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We present a detailed single crystal x-ray difltraction study of twinned orthorhombic perovskites. 
The diffraction pattern can be indexed on a 2apx2apx2ap lattice, but does not obey cubic symmetry 
relations. The data can be modelled in space group Prima with a twin based on a distribution of 
the b axis over three perpendicular cubic axes. This model allows full structure determination in 
the presence of up to six twin fractions from single crystal x-ray diffraction data. 



I. INTRODUCTION 

Transition metal oxides with the perovskite structure 
exhibit a large variety of interesting physical phenom- 
ena, including high Tc superconductivity, ferroelectricity, 
colossal magnetoresistance, and a variety of spin, charge, 
and orbital orderings. Doping can change the magnetic 
and electronic properties, which can be reflected in the 
details of the crystal structure. We are particularly in- 
terested in the structural response of phenomena of elec- 
tronic origin, like Jahn- Teller distortions or orbital order. 

Since the metal-oxygen bond lengths and angles deter- 
mine the exchange interactions, a thorough knowledge 
of the structure is necessary to understand the physics 
behind these phenomena. Typically crystal structure 
information is generated by neutron powder diffraction 
(NPD). NPD gives direct information about the distances 
of lattice planes and is very accurate in determining the 
changes of the lattice parameters with temperature and 
pressure. Using a full pattern Rietveld Analysis the 

atomic positions can be determined with an accuracy of 
10-^-10^3 

The alternative is single crystal X-ray diffraction 
(SXD). SXD is a very powerful experimental application 
to study the crystal structure of a large variety of ma- 
terials. In contrast with powder diffraction, which maps 
all crystal planes to one dimension (d spacing, or angle 
29), SXD observes all crystal plane reflections separately, 
even if they are related by symmetry. The relative co- 
ordinates of the atoms in the unit cell can be determined 
from the intensity distribution with the same accuracy as 
NPD. The major difference is that NPD is sensitive to the 
nucleus, while SXD measures the electron density, which 
is of interest to understand the electronic properties. 

However, for many perovskites SXD is not used be- 
cause of unconventional twinning of the crystals. Twin- 
ning can complicate the structure determination by cre- 
ating extra reflection spots and/or superimposing reflec- 
tions on top of each other. Twinning of the perovskites 
is complex as the cubic parent structure allows not only 
simple twin axis and planes, but also 3D twin relations. 
In this paper we will show a method to analyse twinned 
single crystals. Moreover, this twin model allows a full 
structure determination including detailed information 



on distortions of both structural and electronic origin. 
Besides, we also show that one does not need a larger 
data set, as is usual in the case of twinning. 

Using SXD, we observe a 2apx2apx2ap unit cell. It is 
commonly accepted that the unit cell is •\/2apx2apX-\/2ap. 
The observed unit cell originates from a three- 
dimensional type of twinning that is not restricted to 
manganites. The model is most likely of general appli- 
cation for a large variety of perovskite Pnma crystals. 
This twinning is unique and it involves a distribution of 
the b axis over three perpendicular cubic axes. 

II. THE TWIN STRUCTURE FOR MANGANITES 

The manganites have generated considerable attention 
because of the colossal magnetoresistance effect. The role 
of magnetic order has been widely discussed as the double 
exchange, inducing the ferromagnetic order, is required 
to generate a metallic ground state. The role of orbital 
order is much less understood. While local Jahn- Teller 
distortions are crucial in explaining the localization of 
the charge carriers in the paramagnetic state, the long 
range Jahn- Teller ordering is not well studied, except for 
undoped LaMnOa. We will show that the destruction of 
long range orbital order is a second prerequisite for the 
metallic ground state. In fact, we will show in a separate 
paper that with increasing hole doping of LaMnOa not 
the exceeding of a critical concentration but the suppres- 
sion of the orbital order generates the metallic ground 
state. The orbital ordering in perovskites with degen- 
erate Cg electrons can be easily measured, whereas for 
degenerate t2g electrons the Jahn- Teller distortions are 
much smaller. Furthermore, the Pnma symmetry can ac- 
commodate not only the Jahn- Teller ordering, but also a 
3D rotation of the octahedra, known as the GdFeOa dis- 
tortion. In this paper we will focus on the main reason 
why SXD has not been widely used for these perovskites, 
namely twinning. Twinning in doped LaMnOa originates 
from the transition of the highly symmetric cubic parent 
structure to the orthorhombic symmetry, that accommo- 
dates both the GdFeOa and the Jahn- Teller distortion. 
In thin films ac twinning has been observed. The solution 
of the twin relations in the crystals allow us to study in 
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detail the ordering of these compounds as influenced by 
temperature, magnetic state and doping concentration. 
In this paper we will focus on the twin relations. 

At high temperatures most AMnOa are cubic with a 
unit cell of a = 3.9 A. Due to both the small average 
A-site radius and the JT effect the MnOg octahedra are 
rotated and distorted at lower temperatures. We will 
show that our crystals are twinned. We present here an 
accurate description of the origin of the twinning, the 
detwinning process and show that we can still determine 
the temperature dependence of the structure by measur- 
ing only the reflections of the main fraction in one octant 
of the hkl space. 

Twinning is reported before in AMnOa perovskites 
and related structures. For instance, neutron powder 
experiments on LaMnOs showed reflections, at temper- 
atures above the Jahn- Teller transition, that could be 
indexed on a double cubic unit cell. But full pattern 
refinement was only possible in the orthorhombic space 
group Pnman Single-crystal electron diffraction showed 
a double cubic uuit cell for SrSnOa, even though CaSnOs 
is orthorhombic.B Electron microscopy, however, showed 
the existence of coherent domains, with perpendicular 
orientations of the doubled b axis. The structure was 
deduced via the O'Keefe-Hyde relations J3 

The single crystal was mounted on an ENRAF- 
NONIUS CAD4 single crystal diffractometer. The tem- 
perature of the crystal was controlled by, heating, a con- 
stant nitrogen flow. Initial measurements were done at 
180 K. Temperature dependent measurements were per- 
formed on a diffractometer with an adjustable tempera- 
ture set-up, between 130 K and 300 K. 

III. CRYSTAL STRUCTURE 

Most structure research of perovskites focuses on the 
Mn-0 distances and the Mn-O-Mn angles, as these pa- 
rameters determine the super- and double exchange in- 
teractions. Here, we like to stress the importance of a 
complete structure refinement, including the La-position 
and the rotation/distortion of the MnOg octahedra. The 
basic deformations that determine the deviation from the 
cubic structure are the GdFeOs distortion and the Jahn- 
Teller distortion. 

The basic building block of the perovskite structure is 
a 3.9 A cube with Mn in the centre and O at the face- 
centers. The oxygen ions co-ordinate the Mn to form 
MnOg octahedra. The A atoms are located at the cor- 
ners of the cube. The undistorted 'parent' cubic struc- 
ture rarely exists, but distorts conventionally to an or- 
thorhombic or rhombohedral symmetry. For an ideal per- 
ovskite the ratio between the radii of the A site ion and 
the transition metal ion is such that the tolerance factor 

V2((rM„4/3+) +ro2-) 



is equal to one. In the La-Ca system the tolerance fac- 
tor varies from 0.943 to 0.903 going from LaMnOa to 
CaMnOa. 

A. GdFeOs rotation 

Due to the small radius of the A-site ion, with respect 
to its surrounding cage, the MnOg octahedra tilt and 
buckle to accommodate the lanthanide. This is known 
as the GdFeOs distortion. The cubic state allows one 
unique oxygen position. Due to the GdFeOs distortion 
we need two inequivalent positions to describe the struc- 
ture. 01 is the in-plane oxygen, on a general position, 
X, y, z. Two opposite Mn-01 bonds have the same length, 
but the perpendicular bonds need not to be equal. 02 
is the apical oxygen, located on a fourfold x, j,z posi- 
tion on the mirror plane. Mn-02 bonds are always of the 
same length. Both in the undistorted and the distorted 
perovskite, the 0-Mn-O bond angles are 180° (90°), but 
due to the buckling Mn-O-Mn bond angles are no longer 
180°. A pure GdFeOs distortion can be obtained with 
equal Mn-0 bond lengths. 

B. Jahn- Teller distortion 

The Jahn- Teller effect originates from the degenerate 
Mn ion. Two, possible, distortions are associated 
with the Jahn- Teller effect. Q2 is a orthorhombic dis- 
tortion, with the in-plane bonds differentiating in a long 
and a short one. Q3 is the tetragonal distortion with the 
in-plane bond lengths shortepicig and the out-of-plane 
bonds extending, or vice versaoB. The main result of the 
JT effect, Q2 and Q3, is that the Mn-0 distances be- 
come different, which lifts the degeneracy of the t2g and 
eg levels. 

C. Glazer's view on the octahedra 

There are three possible rotations for a rigid MnOg oc- 
tahedron. One, riy, changes the atomic positions of the 
02 atom only. The others have an effect on the atomic 
positions of both 01 and 02. These movements will also 
affect the position of the La atomQ. GlazeiQ identified all 
possible sequences of rotations for the perovskite system. 
Two of these sequences create the symmetry elements 
that make up the standard perovskites. In Glazer's nota- 
tion: a^b^a^ ^ yields orthorhombic Pnma and a~a+a^ 



*^ p~q'^r'^ means that along [100] a rotation of size p is 
alternated positive and negative, along [010] a rotation of size 
q is always (arbitrarily) positive, and along [001] there is no 
rotation. 
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will result in rhombohedral i?3c. As Glazer worked in 
a pseudocubic 2apx2apx2ap system we have to translate 
these rotations to the orthorhombic unit cell, which gives 
a~6+0. This means that the rotation around the c axis, 
rjz, is zero. This is in agreement with the observation that 
the four positions of 02 near the Mn at 0,0,0 have the 
same distance to the y — plane, albeit two are positive 
and two are negative. By rotating the MnOg octahedra 
around the x axis, two oxygen are raised from their cubic 
position and two are lowered. Here we neglect the influ- 
ence of a rotation rjy on these coordinates. As a result 
oi rjz = the fractional x co-ordinate of 01 is also zero. 
In practice small deviations are found, indicative of the 
non-rigid behaviour of the MnOe octahedra. Typically 
in AMnOa two rotations and the Q2 JT distortion are 
observed. 

One effect of the GdFeOs distortion is that the num- 
ber of formula units per unit cell is enlarged. The unit 
cell is doubled in the b direction, with respect to the 
original cubic cell. The ac plane is also doubled, re- 
sulting in 6 ft! 2ap, and a « c « ^/2ap. Despite the 
rotations and distortion, the structure remains in origin 
cubic. The shifts of the atoms in the unit cell are small. 
Furthermore, Mn stays roughly octahedrally surrounded 
although the Mn-0 distances tend to differ and some O- 
Mn-0 angles may be no longer perfectly 90° or 180°. The 
small deviations from the cubic symmetry are reflected in 
the intensities of the reflections. To illustrate this point, 
we transform an arbitrary, cubic, reflection hkl to the 
orthorhombic setting by h'k'l' = h+l,2k,h-l. Thus or- 
thorhombic reflections with k' even and h' and /' both 
even or both odd stem from planes that already existed 
in the original cubic unit cell. We observe in diffraction 
patterns of orthorhombic perovskites that reflections that 
originate from cubic crystal planes have higher intensities 
than those that do not. 



IV. TWIN MODEL 

Reflections in reciprocal space are experimentally ob- 
served at a regular distance in three orthogonal direc- 
tions, corresponding to a cubic lattice spacing of 7.8 A 
in real space. Although the three orthonormal axes have 
equal lengths, we could not observe the threefold rota- 
tion axis along (111), required for cubic symmetry. How- 
ever, studying the intensity distribution of planes in hkl 
space with constant h, k or I showed much regularity 
as is shown in Fig. |l|. We propose a twinned struc- 
ture consisting of coherent Pnma domains. Due to par- 
tial overlap this results in a metric cubic system with 
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FIG. 1. The number of reflections of the different types 
versus intensity ranges. Low intensities, I < 200 are mostly 
EOO, while high intensities, I > 2500, are of the types eee 
and EEO. 

Twinning is often observed in crystals with a reduced 
symmetry. For instance, the orthorhombic perovskite 
LaMnOs has lower symmetry than cubic SrTiOa. Con- 
ventional twin models keep one characteristic axis un- 
changed and form the domains by rotation around that 
axis. This is commonly observed in constrained epitax- 
ial thin films and pseudo-two-dimensional crystals like 
YBa2Cu307_5. Another standard twin is inversion twin- 
ning, found in non-centrosymmetric systems. Due to the 
inversion twin, pscudo-centrosymmetry is obtained or a 
net polarisation in ferroelectric compounds is reduced 
to zero-u'EI Both twin models have usually the property 
that either all twin domains have reflections that lie on 
top of each other, merohedral twins, or there is a non- 
commensurate set of reflections, e.g. in monoclinic unit 
cells. Here we propose a more extensive form of twin- 
ning, where the different reflections of different domains 
partially coincide, giving rise to the observed metrically 
cubic system. 

The transformation of cubic to orthorhombic symme- 
try requires a designation of a, b, and c with respect to 
the degenerate cubic axes. There are three possibilities to 
position the doubled b axis along the three original cubic 
axes. Thus we propose that the three fractions' b axes are 
oriented along the three original axes of the cubic unit 
cell, as sketched in Fig. ^. This twin model consists of a 
Pnma unit cell, transformed by rotation along the 'cu- 
bic' [111] axis. We still have the freedom to choose the 
a and c axes, perpendicular to the b axis, and rotated 
45° with respect to the cubic axes. Therefore, this model 
yields six different orientations of the orthorhombic unit 
cell. As the differences between a/^/2, b/2 and c/V2 are 
small, we observe the reciprocal superposition of the six 
orientations as metrically cubic. 
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FIG. 2. Two of three possible orientations of the Pnma 
unit cell in the observed 2apx2apx2ap {a"xb"xc") unit cell. 
b'l , b'2 and 63 denote the orientations of the doubled axes along 
the original cubic axes, ap. Note that the choice of a' and c' 
is still open, after choosing an orientation for the b' axis. 

The standard refining programs can work with twin 
models, but only if they consists of merohedral twins. In 
our case, the reflections do not correspond with an or- 
thorhombic unit cell, but with a, twice as large, cubic 
unit cell. In appendices |^ and we will elaborately 
explain the transformation from orthorhombic to double 
cubic and the way we worked around the refining pro- 
gram. 

We first measured a hemisphere in hkl space with -20 < 
h < 20, -20 < A: < 20, and < Z < 20. On this dataset, 
we could refine our twin model with the six fractions, 
including the Ca concentration on the A site. Measuring 
a large range in hkl space requires a large amount of time, 
roughly 10 days in the present set-up. Conventionaly, 
for crystals with orthorhombic unit cells measuring one 
octant is sufficient. From refinements on selected parts of 
the dataset, we concluded that we could investigate the 



structure at different temperatures by only measuring the 
positive octant of the main fraction, thereby limiting the 
measuring time considerably. 



V. DISCUSSION 

The observed refiections could be indexed on a cubic 
lattice with a ~ 7.8 A. Although the three axes had equal 
lengths, we could not observe the required threefold sym- 
metry axis along (111) for cubic symmetry. This is shown 
for EEO reflections in Fig. ||. Furthermore, studying the 
intensity distribution of hkl planes with constant h, k 
or I showed much regularity, as is shown in appendix 
0. Rodrigruez-Carvajal et al. noted that their neutron 
powder spectra of pure LaMnOs could be indexed in a 
cubic 2apx2apx2ap mat cell, above the Jahn- Teller tran- 
sition temperature.cJ They did not observe a splitting of 
the peaks. Nevertheless, they could not refine the struc- 
ture in a cubic unit cell. Note that in powder diffraction 
spectra there is no direct method to observe the three- 
fold symmetry along the body diagonal. The refinement 
of the spectra has been done in the conventional Pnma 
setting. They suggested that the powder is twinned but 
could not give evidence for that, as they only studied 
powders. The observed reflections, the systematic ex- 
tinctions and the non-cubic intensity distribution suggest 
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FIG. 3. The number of reflections versus the types OEE, 
EEO and EOE for different intensity ranges. OEE, b' paral- 
lel to a", has mostly very low intensity. Reflections EOE, b' 
parallel to b" , have a broad range of intensities. In cubic sym- 
metry, these distributions should be identical for absorption 
corrected data, which we considered. 

Typically, Ca-doped LaMnOa has space group Pnma. 
We propose a twinned structure consisting of coherent 
Pnma domains, which results in a metric cubic system 
with a « 7.8 A. The model is identical to the. twin- 
ning observed by electron microscopy in CaSnOs.H They 
could however not refine their 2apx2apx2ap data on dif- 
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ferent Pnma domains, ancLdeduced the structure from 
phenomenological relationsa and the similarities with 
SrSnOs- In the present work, the different contributions 
of the twin fractions to the total integrated peak intensity 
are taken into account. The refined model consisted of 
the regular parameters in single crystal diffractometer. 
We considered for each observed reflection the related 
orthorhombic reflections for the applicable twin fractions 
and used an identical unit cell. This way, we can refine 
both the atomic positions of the asymmetric unit and the 
volume ratio of the twin domains simultaneously. 

Several crystals were measured and refined with this 
detwin model. The refinements showed that for every 
crystal the distribution of the volume over the twin frac- 
tions is different. This is an indication that our crystals 
are really twinned. Other structural deformations would 
lead to a constant, not sample dependent, effect on the 
structure factors and therefore to the same fractions of 
the twin domains. We can also conclude that the size of 
the twin domains must be slightly smaller than the mag- 
nitude of the measured crystals, i.e. several tenths of a 
millimeter. Larger twin domains would give rise to crys- 
tals of one single domain and smaller domains are more 
likely to produce a constant spreading of the domains. 
Non-regular, though constant, distributions could signal 
a preferential growth direction of the crystal. 

A. Intensity distribution 

Careful analysis of the intensity distribution showed 
extra evidence for the proposed models. The full dataset 
shows that reflections hkl of type eee^ arc by far the 
strongest. Reflections having one odd Miller index, OEE, 
EOE, EEO, are second in intensity, eoe has the highest 
intensity of these reflections as seen in Fig. |. To under- 
stand the intensity distribution in the 2apx2apx2ap unit 
cell, we examine the structural deformations. Despite 
the effect of Jahn- Teller distortions, GdFeOs rotations 
and twinning, the structure remains in origin cubic. This 
means that the intensity distribution of the peaks with 
high intensity will always mimic the intensity distribu- 
tion of the peaks of the undistorted cubic perovskite. An 
arbitrary reflection, in the cubic unit cell, hkl, is trans- 
formed to the orthorhombic unit cell as h + I, 2k, h — I. 
Thus orthorhombic reflections with h' and both even 
or both odd and k' even stem from the original cubic unit 
cell. Although the fact that the h'k'l' =eoe reflections 
do not originate from the original cubic planes, they are 
allowed in the Pnma symmetry. Reflections h'k'l' that 
satisfy these conditions have higher intensities than those 
that do not. 



^'E denotes an even value for the Miller indices, O stands for 
odd 



To double the Pnma unit cell to a 2apx2apx2ap unit 
cell, a reflection h'k'l' is transformed to h"k"l" as h' + 1', 
k', h' — I' . Please note that we do not explicitly mention 
the different orientations of the orthorhombic unit cell. 
If the regarded b axis is not parallel with k" , than the in- 
dices should be cycled appropriately. The transformation 
yields that double cubic reflections with k" = k' = 2n 
originate from single cubic reflections, with k — n. How- 
ever, k" = 2n + 1 reflections originate from single cubic 
non-integer indices k, and will therefore have less inten- 
sity. This explains why the measured h"k"l" =eee reflec- 
tions are strongest and eoe are one order of magnitude 
less strong. We have seen that h" = h' + l' and I" — h' — l' 
thus h" and I" are always both even or both odd. There- 
fore EOE has only contribution of the fractions with k" 
corresponding to the orthorhombic b axis. OEE and EEO 
originate from fractions with their b axis parallel to h" 
and I", respectively. The fact that our observed eoe is 
stronger than OEE and eeo, indicates that, by chance, 
the observed b axis, or k" direction, is parallel with the 
b axis or k' direction of the largest twin fraction. 

B. Refinement 

We determined both the Ca content and the volume 
fractions of the twins by measuring a full hemisphere of 
hkl space. We compared the results with a reflnement 
on only one octant of these measured reflections. The 
outcome was equivalent within the error bars. We re- 
fined again fixing the twin fractions and Ca concentration 
to the values found for the first refinement on the large 
dataset, and allowing only the atomic co-ordinates and 
anisotropic displacement parameters (adp's) to change. 
This resulted in the same structural model as found 
with the complete structure determination on the largest 
dataset. We concluded that apparently the refinement 
with fixed fractions and Ca concentration is insensitive 
for the size of the dataset. This can be understood if 
we view the symmetry relations and Friedel pairs. Due 
to the Patterson symmetry, mmm the relation F{hkl) = 
F{hkl) is valid for all structure factors. Friedel pairs are 
refiections that have intrinsically the same structure fac- 
tor, related thus by symmetry and ignoring absorption 
and anomalous scattering effects. Although we have six 
different orientations, one octant chosen for a particular 
orientation will still contain data of all fractions that can 
be transformed in such a way that every h'k'l' is included 
for all three b directions. 



VI. CONCLUSION 

We have shown that crystals of Lai-ajCa^^MnOa, with 
X — 0.19, have the orthorhombic Pnma space group, al- 
though they appear in single crystal diffractometry met- 
rically cubic. The metric cubic appearance results from 



5 



twinning in which the doubled b axis of the Pnma cell is 
oriented along all three cubic axes of the single cubic unit 
cell during crystal growth or cooling. This twin model 
allows successful structure refinements of single crystal 
x-ray diffraction data. 

APPENDIX A: DATA ANALYSIS 

Here we describe the effect of the twin model on the 
visibility of the standard Pnma reflection conditions as 
well as a detailed analysis of the intensity distribution of 
a model crystal with Pnma symmetry and of the mea- 
sured reflections in double cubic setting. The first step 
in space group determination is to look for systematic 
absences in the list of reflections. We observe no inten- 
sity for the reflections shown in table |. Note that the 
indices are in double cubic setting and therefore allow 
cyclic permutation, i.e. hOO, OfcO and 001 are all rep- 
resented by hOO. These reflection conditions allow only 
the space groups P2i3 and P4232. However, these space 
groups have as reflection condition only hOO : h = 2n. 
This suggests that we have more information than can be 
attributed to these space groups. Attempts to determine 
the structure using these space groups were unsuccess- 
ful. Furthermore, the observed and, for cubic systems, 
anomalous, reflection conditions suggested that we might 
be looking at a twinned crystal. In the next paragraph, 
we present the transformation of the reflection conditions 
for orthorhombic Pnma according to the presented twin 
model. 

TABLE L Extinct reflections as observed in the 
2apx2apx2ap data set. 



Reflection Extinction condition 

hOO hj^2n 

hhO /i / 2n 

hkh h 7^ 2n, k = 2n 

hkh h 7^ 2n, k = 2n 



1. Extinctions 

The Pnma reflection condition Ofc'O : k' — 2n is trans- 
formed to double cubic as Ok"0 : k" = 2n and cychc 
permutation. This corresponds to the reflection condi- 
tion h"00 : h" = 2n in tabic ||. Pnma reflections of the 
type h'Oh' are transformed to h"00, with h" = 2h' . This 
implies that h/Oh' never contributes intensity on reflec- 
tions h"OQ, with h" odd. Conversely, if we do not observe 
intensity on all reflections h"00 : h" ^ 2n, then all of the 
constituting Pnma reflections should be absent. 

To transform and understand the other reflection con- 
ditions of Pnma we have to take into account that an or- 
thorhombic, and also a cubic, unit cell implies that hkl is 
equivalent upon sign reversal of each of the Miller indices. 
The Pnma reflection condition h'OO : h' = 2n becomes 
h"Oh" : h" = 2n and similarly OOV transforms _to V'OF. 
But as we just stated, h"Oh" is identical to h"Oh" . If we 
combine reflections then a reflection condition is fulfilled 
if any of the contributing parts has intensity. To obtain 
extinction, all contributing reflections should be extinct. 
As OOZ' has no extinction condition, i.e. intensity for all 
V, they will contribute to all h"Oh"/h"(W' reflections, 
and the Pnma reflection condition h'OO : h' = 2n is 
masked. 

The overlap of different reflections occurs also for Ok'h' 
and h'k'O. These transform to h"k"h/', and h"k"h". 
Therefore, we cannot disentangle Ok'h' and h'k'O in the 
measured double cubic cell. The Pnma reflection con- 
ditions for Ok'l' and h'k'O are k'+l' = 2n and h' = 2n, 
respectively. h"k"h" and h"k"h" are found to be extinct 
for h" 2n and k" — 2n in the double cubic setting. 
To have these reflection extinct, all contributing Pnma 
reflection conditions must be extinct. This implies that 
both h' = 2n and k' -\-h' = 2n, the reflection conditions 
for the contributing Ok'h' and h'k'O, may not be fulflUcd. 
Thus h" 7^ 2n and this yields k" = 2n. 

Now we consider reflections of the type h"k"h" that 
are not extinct. Our twin model yields six fractions that 
can contribute to the considered reflections. 

1. h" = 2n and k" ~ 2n. Reflections eee have con- 
tributions from the four fractions with the doubled 
6'-axis parallel to either a" or c", this is also true 
for reflections OOO. In addition, we also have con- 
tributions from the two fractions with b' parallel to 
6", i.e. from h'k'O and Ok'h'. 



2. h" ^ 2n and k" = 2n. Reflections OOO, in addition 
to the four fractions with the doubled 6'-axis paral- 
lel to either a" or c", have one additional contribu- 
tion from Ok'h' . In the special case that h" ~ k", 
for both EEE and OOO, then two of the four frac- 
tions, parallel to a and c contribute. 

3. h" = 2n and k" ^ 2n. Reflection eoe, have only a 
contribution from the h'k'O reflection, with b' par- 
allel to the observed b" axis. 
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We have shown that the intensity of observed h"k"h" 
reflections depend on odd or even. We conclude that 
h"k"h" reflections are extinct for OEO, while the other 
reflections have contributions from up to six twin frac- 
tions. 



2. Intensity in double cubic setting compared with 

Pnma 

We can also learn something about the twinned ori- 
gin by studying the intensity distribution. Our twin 
model provides a natural explanation for the hierarchy 
in intensities of the reflections. First, we consider the 
known transition from the single cubic unit cell to the or- 
thorhombic one to elucidate the patterns that are inher- 
ent to perovskites. Then we proceed to show how these 
patterns are convoluted in the double cubic, twinned 
model. 

The intensity distribution of a normal Pnma per- 
ovskite is shown in Fig. H. In conventional orthorhombic 
Pnma the reflections h'Pl' with h' + 1' = 2n are stronger 
than h' + 1' ^ 2n. From these, the k — 2n reflections are 
stronger than the k ^ 2n ones. This can be attributed 
to their origin in the 3.9 A cubic structure. The h'k'l' 
reflections in orthorhombic setting originate from hkl in 
the single cubic cell for h' = h+l and I' = l-h. Therefore, 
h' and V will always be both even or both odd. In other 
words if h'-l' ^ 2n than the reflections originate from a 
crystal plane in the single cubic setting with a non-integer 
Miller index, i.e. a superlattice reflection and their in- 
tensities are usually weak. The same argument can be 
applied to k. As k' = 2fc, thus k' = 2n originate from a 
regular crystal plane, k' ^ 2n is the superlattice reflec- 
tion. These k' ^ 2n however appear to have somewhat 
more intensity than the h'-V ^ 2n reflections. 



100 



10 



CD 



Tk 1 




^ 0.1 
0.01 

FIG. 4. The calculated intensities for a conventional per- 
ovskite, subdivided with respect to odd and even of the Miller 
indices, k = 2n reflections are strongest, if h' + I' = 2n. Re- 
flections with h'+l' 7^ 2n have roughly 100 times less intensity. 

We observe the following patterns in the intensity dis- 
tribution in the double cubic cell. First, reflections eee 
are generally the strongest, EOO usually absent, or very 



weak. If we ignore the allowed permutations in cubic set- 
ting, the following hierarchy can be made, eoe reflections 
were much stronger than eeo, which are again stronger 
than OEE. Likewise, OEO reflections were, although weak, 
stronger than OOE and again stronger than EOO. 

In conventional Pnma structure we found that if 
h'+l' 2n than reflections with k' ^ 2n are stronger 
than those with k' = 2n. In our notation, for the Pnma 
structure this yields (eee=OEO) >> (ooo=eoe) >> 
(00E=E00) > (oee=eeo). So there are roughly four 
intensity groups. We can correlate these four groups 
with the four types of reflections we measured in the su- 
percubic setting, with the following intensity hierarchy: 
eee>>eeo>>ooo>eoo. 

We consider here the transformation from the or- 
thorhombic to the double cubic setting for the in-plane 
indices, h' and I'. Double cubic indices are calculated 
with h'+l' and h'-l', therefore h" and I" will be both 
even or both odd. Reflections of the type E?0 or 0?E will 
have no intensity of the twin with b' parallel to b" , as they 
are not originating from an 'integer' crystal plane. But 
we have a mixing of different orientations of the b' axis 
along the three cubic axes a", b" and c". We can cycle 
the indices of these reflections such that we get either 
E?E or 0?0, and from the corresponding twin, b' parallel 
to a" or c" we do have intensity on these reflections, eoe 
and OEO will only have contributions from this particular 
setting. EEE and OOO can be cycled and in general have 
contributions from all six orientations. 



3. distribution of the twins 

If we consider that the main part of the intensity in 
the double cubic setting originates from the largest twin 
fraction, we can search for the orientation of the largest 
fraction. That b' axis is expected to be oriented along 
one of the measured double cubic axis. We can not dif- 
ferentiate the two fractions with parallel b' and perpen- 
dicular a' and c' . We considered groups of reflections 
with cycled indices, having either one odd, OEE, or one 
even index, EOO. Roughly, the intensities of the reflec- 
tions within these groups occurred with ratio 70:20:10. 
This suggested that one twin had 70% of the volume, the 
others 20% and 10%. Now we sorted these reflections to 

flnd the corresponding orientations. The sort parame- 
ter is the observed double cubic axis that corresponded 
with the odd indices in OEE, and vice versa. Only cy- 
cled reflections that occurred three times with measur- 
able intensity were taking into account. We found that 
la" ■ h" ■ Ic" = 5 : 80 : 15. The constraint that all 
three intensities have positive values, ignores the weak- 
est reflections. It also ignores Bijvoet pairs, as we only 
measured —20<l"<2. If we indeed had a cubic system 
then these variations should be zero within standard de- 
viation. The distribution of the intensities with respect 
to the different axis, strongly suggests that the crystal is 
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twinned. 



APPENDIX B: IMPLEMENTATION IN 
REFINEMENT 

Here we describe the model and the appUcation by 
SHELXL. The standard refining programs can work with 
merohedral twin models. However, in our case, the re- 
flections do not coincide with one orthorhombic unit cell. 
Every observed reflection has contributions of up to six 
twin fractions. We measured both in the double cubic 
setting as in the orthorhombic setting of the main twin 
fraction. We transformed the h"k"l" indices to the six 
possible twin orientations. Three of the possibilities are 
given by 

h'k'l' = + l"),k", - I") (Bl) 

h'k'l' = + k"),h'\ i(r - k") (B2) 

h'k'l' = ]^{k" + h"), I", i(fc" - h") (B3) 

The other three can be acquired by changing h' for I' and 
V for -h'. 

A new software program TWINSXL was developed to 
transform the standard data file, HKLS, by using the 
appropriate transformation matrices from a second in- 
put filellj The new data file constitutes of lines of ft-, k, I, 
intensity plus standard deviation and twin fraction num- 
ber. We used the "HKLF 5" option of SHELXL to reflne 
data. The refinement uses a crystal model for the or- 
thorhombic structure with the normal, adjustable vari- 
ables. Five variables for the twin fractions were added, 
the sixth fraction is calculated as the complement of the 
other five fractions. The sum of the appropriate calcu- 
lated intensities for all fractions was compared with the 
observed integrated intensities. 
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